function [ s ] = boundaryQuadrature( p0, q0, q1, qn )
%BOUNDARYQUADRATURE Summary of this function goes here
%   Detailed explanation goes here

d = q1 - q0;

l = sqrt(d(1)^2 + d(2)^2);

v1 = evaluateQuadraturePoint(p0, q0, qn);
v2 = evaluateQuadraturePoint(p0, q1, qn);

s = (v1 + v2) * l / 2.0;

end

